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We study the effects of field modulation on the energy spectrum of an electron in a two-dimensional 
bipartite periodic lattice subject to a magnetic field. Dependence of the energy spectrum on both 
the period and the strength of field modulation is discussed in detail. Our main finding is that 
introducing field modulation drastically changes the energy spectrum and the localization properties 
of the system appearing in the absence of field modulation; the degeneracies induced by a uniform 
magnetic field are broken and the resultant energy spectrum shows a dispersive band structure, 
indicating that most of Aharonov-Bohm cages become unbounded. The effects of field modulation 
on the superconducting transition temperature and the critical current in a wire network with the 
same geometry are also discussed. 
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I. INTRODUCTION 

The physics of magnetically induced frustration in vari- 
ous two-dimensional (2D) structures including square, 1-3 
rectangular, 4,5 triangular, 6 honeycomb, 7,8 aperiodic, 9 
quasiperiodic, 10,11 fractal, 12 and even random 13 geome- 
tries has attracted much interest in condensed-matter 
physics for several decades. Recently, Vidal and co- 
workers 14 presented a new localization mechanism in- 
duced by a uniform magnetic field for nonintcracting elec- 
trons in a 2D bipartite periodic hexagonal structure (the 
so-called T3 geometry), where the unit cell contains three 
sites, one sixfold coordinated (called the 'hub' site) and 
two threefold coordinated (called the 'rim' sites). Within 
the tight-binding (TB) approximation, they showed that, 
due to fully destructive quantum interference, the eigen- 
states at half a magnetic flux quantum per elementary 
rhombus (briefly, half a flux) are extremely localized and 
bounded in Aharonov-Bohm (AB) cages. Very recently, 
Abilio and co-workers 15 performed transport measure- 
ments on a superconducting wire network with the T3 
geometry and confirmed the field-induced localization ef- 
fect by observing a depression of the superconducting 
transition temperature T c and the critical current J c at 
half a flux. 

Since the calculations in Ref. 14 were performed un- 
der the basic assumptions of (i) the perfectness and the 
infinite size of the T3 geometry and (ii) the spatial uni- 
formity of the magnetic field, it may be very interesting 
to address a question of what happens under the situ- 
ations beyond the two assumptions. Of course, effects 
beyond the assumption (i) were already discussed by the 
authors of Ref. 14; the randomness (such as random mod- 
ulations of hopping terms and fluctuations in the tiling 
areas and in the transmission matrix along the edges), 
which is inevitable in real systems, was expected to al- 
ter or even destroy the phase matching essential for the 
field-induced localization effect. In addition, the authors 
of Ref. 15 argued that the incomplete suppression of the 
experimentally observed J c at half a flux might be at- 



tributed to the network's finite size. However, effects 
beyond the assumption (ii) are not examined yet and re- 
main an open question. Thus, in this paper, we would 
like to address the question and clarify the effects of a 
nonuniform magnetic field on the energy spectrum and 
the transport properties of the T3 geometry at rational 
fluxes, especially at half a flux. As a concrete and sim- 
ple example, we consider a magnetic field with a periodic 
modulation and investigate the effect of field modulation 
on the stability of AB cages and on the characteristics of 
T c and J c . 

The contents of this paper are organized as follows. An 
eigenvalue equation taking into account field modulation 
is derived in Sec. II. Analytic and numerical results for 
the effects of field modulation on the energy spectrum 
at rational fluxes and the localization properties of the 
eigenstates at half a flux are presented in Sec. III. The 
effects of field modulation on T c and J c of a supercon- 
ducting wire network are also discussed in this section. 
Finally, Sec. IV is devoted to a summary. 



II. FIELD MODULATION AND THE 
EIGENVALUE EQUATION 

We consider an electron in a 2D rhombus tiling under 
a spatially modulated magnetic field as 



B=[B + B mod (x)} z, 



(1) 



where Bq (B moc i) denotes the uniform (modulated) 
part of the applied magnetic field. Among possible 
types of modulated fields, we pay attention to the one- 
dimensional (ID) sine-modulated field as 



B m od(x) = Bi sin 



27ra; 



(2) 



where T x is the period of the modulation along the x 
direction. Under the Landau gauge, the vector potential 
is given by 
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The TB equation that describes an electron on a 2D 
lattice subject to a magnetic field reads 



(4) 



where tij is the hopping integral between the nearest- 
neighboring sites i and j, and the phase factor 7^ is 
given by 
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4>a = hc/e being the magnetic flux quantum. For the 
sake of simplicity, we set Uj = 1 . The phase factor along 
the x direction in Fig. 1 is zero under the Landau gauge. 
Denoting the phase factor along the upward direction on 
the line 1(2, 3, 4) in Fig. 1 as 71(2,3,4), and using the wave 
function at a hub site as ip{x 7 y) — e' lkvV ip(x), Eq. can 
be written as 
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where n = ^fik y a/2. 

Denoting ip m — ip{x) at x — 3ma/2, the phase factors 
can be written as 
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4>{= V3B a 2 /2) being the uniform background magnetic 
flux through the elementary rhombus. Now let us define 
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Then, Eq. (||) can be reduced to a ID equation, which is 
a generalized version of the eigenvalue equation derived 
in Ref. 14: 



\lp m = A m+1 tfj m+1 + C m 1p m + A m t/j m -l, 



where 



A m = cos(> m ) cos(^+ + k), 
C m = cos(v m ) cos(i^+ + 2k), 
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A close inspection of Eqs. (|llj) - ( |13| ) for a reduced 
rational flux / = p/q with mutual primes p and q enables 
us to write the Bloch condition along the x direction as 



Here, 77 = 3k x a/2 and N is given by 
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where 
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q' ', p' ', and being integers. Using Eqs. ( |ll| ) an( l fll4|), 
we obtain the characteristic matrix as 
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Thus, by diagonalizing Eq. (|17|), we can obtain 2N energy 
eigenvalues {±(6 + 4Ai) 1 / 2 ; i = 1,2,..., N} for a given fc 
and the full energy spectrum by sweeping all the k points 
in the magnetic Brillouin zone (MBZ). 



and 



A = 



E 2 -6) 



(10) 



III. RESULTS AND DISCUSSION 
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A. E — f diagram 

When a periodic field modulation is introduced, a 
new length (i.e., T x ) is added into two characteristic 
lengths, the lattice constant a and the magnetic length 
I = (H/eBo) 1 / 2 . Thus, a subtle interplay among them 
can change the energy band structure and the localiza- 
tion properties of the system obtained in the absence 
of field modulation. We plot in Fig. 2 the E — f dia- 
grams without and with field modulation. In the calcu- 
lations, we took T x = 3 and q = 101 (1 < p < 100), 
and swept (20 x 20) k points in the MBZ. Comparing 
Fig. 2(b) with Fig. 2(a), we observe the following effects 
of field modulation. First, the occurrence of overlapping 
between neighboring subbands makes most of small gaps 
closed. This phenomenon of subband broadening (or gap 
closing) clearly appears for large values of (3 and/or /. 
Second, introducing field modulation lowers the symme- 
try of the energy spectrum. The translational symmetry 
[i.e., E(f+1) = E(f)] and the reflection invariance about 
/ = 1/2 [i.e., E(f) — E(l — /)] appearing in a uniform 
magnetic field are no longer held in the presence of field 
modulation. Note, however, that the energy spectrum 
still exhibits the symmetry with respect to E = and 
the reflection symmetry about / = 0. The former re- 
flects the bipartite geometry of the rhombus tiling and 
the latter implies that there is no way for an electron to 
discern the direction of the magnetic field. 



B. Energy dispersion at half a flux 

Before proceeding further, it may be worth noting that 
there are two types of localization effects in the T 3 ge- 
ometry. One is the topological localization appearing at 
E = 0, which originates from the local topology of the 
lattice. 16 Since this localization is well understood and 
independent of the magnetic field, we omit further dis- 
cussion on this phenomenon. The other is the uniform 
field-induced localization appearing at E = ±V6 at half 
a flux. 14 We pay attention to how this kind of localization 
is influenced by introducing field modulation. Note also 
that, owing to the reflection symmetry about E = 0, we 
pay attention to the energy spectrum only with E > in 
the discussion below. 



1. Case ofT x =2 

In this case, TV = 4 and direct diagonalization of the 
(4 x 4) Hamiltonian matrix yields the following energy 
dispersion: 



A± = {{A\ + At + A 2 /2j 2 ± [Al(A 2 
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A close inspection of Eq. (|18|) shows that the upper and 
lower edges of the energy dispersion with E > are given 
by 
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and the bandwidth is given by 
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Equation ( p2| ) clearly shows that the bandwidth varies 
from to 4V3 with varying (3. Note that this is a re- 
markable effect of field modulation; the modulated field 
breaks the degeneracy at E = \/6, which appears in a 
uniform magnetic field and makes the energy spectrum 
dispersive. 



2. Case ofT x =3 

By means of a similar method for the case above, we 
have the energy dispersion 
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where 



f(k x , k y )— cos 2 ?ycos 2 n + sin 2 r\ sin 2 n, 

g(k x , k y )= sin 2 r\ cos 2 n + cos 2 77 sin 2 k (25) 



and 
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An inspection of Eq. (|24|) shows that the upper and lower 
edges of the energy dispersion with E > are given by 



E u = [6 + 4v/2~D^ V2 , E d =\&- A^2D~ 2 ] ^ , (27) 
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and the bandwidth is given by the same form as Eq. ( |22| ) 
with replacing D\ by D2, which varies from to 4-\/2 with 
varying (3. We plot in Fig. 3 the energy spectrum as a 
function of (3, where the vertical lines are results obtained 
by numerically diagonalizing Eq. ( |l7j ) and the boundary 
curves are obtained directly from Eq. The inset 

shows the (3 dependence of the energy spectrum up to 
P = 6. The formation of a dispersive energy spectrum 
instead of a highly degenerate point spectrum can be 
clearly seen in the figure. 

Another remarkable effect of field modulation can be 
found in the density of states (DOS). We plot in Fig. 4 the 
DOS for several values of (3. The DOS for (3 = consists 
of a 5-function peak at E — -\/6, as elucidated in Ref. 14. 
However, when /3 ^ 0, the DOS at E = V6 vanishes 
and E = \f& becomes an edge that connects the two 
subbands whose DOS's exhibit the well-known pagoda 
shapes with a logarithmic singularity in the middle of the 
subbands. The formation of a dispersive energy spectrum 
and the vanishing of the DOS at E = \/6 indicate that 
the modulated field makes most of AB cages unbounded, 
which in turn implies that the transport properties of the 
system at half a flux will exhibit quite different behaviors 
from those under a uniform magnetic field. 

3. Case ofT x >4 

As an example to illustrate the fact that the energy 
spectrum sensitively depends not only on (3 but also T x , 
we plot in Fig. 5 the energy spectrum for T x = 4 as a 
function of (3. In this case, N — 8 and there are eight 
subbands with E > 0, some of which may overlap or may 
have finite gaps between them depending on f3. Figure 5 
shows that the number of distinguishable subbands runs 
from 1 to 6 with increasing f3 up to 1. We found that 
the occurrence of subband splitting (or gap opening) is a 
generic feature of the energy spectrum at half a flux for 
T x > 4. We also found that the subband splitting occurs 
more easily when T x ^ 3T X than when T x = 3T X . 

The T x dependence of the DOS also exhibits an inter- 
esting feature. For illustrative purposes, we plot in Fig. 6 
the DOS for T x = 4 for several values of (3. It can be seen 
that E — V6 locates near an edge of a subband for small 
values of f3, while it locates in a subband for large val- 
ues of (3. Note, however, that even in the latter case the 
DOS at E — \/Q is still negligibly small compared with 
the integrated DOS and will have a very little effect on 
the transport properties of the system, if any. 

C. Energy dispersion at generic rational fluxes 

The (3 dependence of the energy spectra for generic 
rationals / exhibits a behavior similar to the case of 
/ = 1/2; the phenomena of subband broadening and gap 
opening are generic features under the modulated field. 



We show an example in Fig. 7, where the energy spec- 
trum with E > is plotted for / = 2/3 and T x = 4. 
When [3 — 0, the energy spectrum consists of two sub- 
bands that touch each other at E = 0, resulting in a 
gapless single-band structure. 1 However, when the field 
modulation is turned on, a gap opens between the two 
subbands. Besides, with increasing (3, there occurs fur- 
ther splitting of each subband into several sub-subbands 
as well as the gap closing and reopening. 

D. Effect of /3 on Tc in a wire network 

Now we are in position to discuss the effects of field 
modulation on T c and J c in a wire network. In the case 
of the T3 geometry, T c (/) is directly related to the edge 
eigenvalue of Eq. ( |Tl| ) by 2,15 

1 -TM =c r C0S {-vm)\ ' (28) 

where C is a constant that is proportional to the square 
of the superconducting coherence length at zero temper- 
ature and E e (f) is the edge eigenvalue at a rational flux 
/. For simplicity, we set C = 1. 

We plot in Fig. 8 the phase boundaries for a system 
with T x = 3. In the calculations, we took q = 101 
(1 < P < 100). It can be clearly seen that the phase 
boundary for (3 = is symmetric about / = 1/2 and 
distinct downward cusps occur at low order rationals 
/ = 1/3,2/3,1/6,5/6,2/9, and 7/9. Besides, an upward 
cusp exists at / = 1/2, which reflects the field- induced 
localization properties of the eigenstates at half a flux, 
as demonstrated in Ref. 15. However, when the field 
modulation is introduced, the phase boundary exhibits 
several distinctive features from the case of (3 = as 
follows. First, the phase boundary becomes asymmet- 
ric with respect to / = 1/2, which reflects the break- 
ing of the reflection invariance of the energy spectrum 
about / = 1/2. Second, the phase boundary for / > 1/3 
severely changes in a nontrivial way even for small values 
of (3, while the boundary for / < 1/3 is less influenced 
even by large values of (3. To put it concretely, of the 
cusps at strong commensurate fields exhibited in the ab- 
sence of field modulation, the cusps at / = 1/3 and 1/6 
are still clearly visible, while the cusp at / = 5/6 becomes 
invisible by introducing field modulation. Meanwhile, 
the cusp at f — 2/3 is visible (invisible) at (3 = 0,0.6 
(f3 = 0.3). Third, with increasing f3, the T c (/) depres- 
sion up to / < 4/5 decreases independently of T x . The (3 
dependence of the T c depression at f — 1/2 is shown in 
Fig. 9, where a monotonic decrease of 1 — T c (/)/T c (0) is 
clearly seen. Meanwhile, the amount of the T c depression 
for / > 4/5 is inconsistent with the modulation strength. 
Fourth, the cusp at / = 1/2 moves downward by intro- 
ducing a field modulation, which is the most remarkable 
effect of field modulation on T c . This phenomenon in- 
dicates that the localization properties of the system at 
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/ = 1/2 becomes similar to those at other rational fluxes 
(such as / = 1/3, 1/2, 1/6, 5/6) in which the eigenstates 
exhibit an extended nature. Fifth, the value of / where 
the maximal T c depression is achieved lowers as (3 in- 
creases. 

Now, let us briefly discuss the effect of field modulation 
on the critical current J c (f), which is closely related to 
the band curvature, d 2 E e (f)/dk 2 , near the band edge. 
When (3 = 0, due to the absence of dispersive states, 
J c (f = 1/2) vanishes completely, as discussed in Rcf. 15. 
However, when (3 ^ 0, the formation of a dispersive band 
structure makes d 2 E e (f)/dk 2 ^ 0, and hence J c (f = 
1/2) will have a finite value, which is another remarkable 
effect of field modulation. 



IV. SUMMARY 

We have studied the effects of field modulation on the 
energy spectrum of an electron in a 2D bipartite peri- 
odic lattice subject to a magnetic field and on the su- 
perconducting transition temperature in a wire network 
with the same geometry. We have shown that the energy 
spectrum sensitively depends on both the period and the 
strength of field modulation. Our main finding is that 
the field-induced localization properties of the lattice at 
half a flux are drastically changed by introducing field 
modulation; the modulated field breaks the degeneracy 
induced by a uniform magnetic field to make the energy 
spectrum dispersive, where the number of distinguish- 
able subbands sensitively depends on both the period 
and the strength of field modulation. The formation of 
a dispersive energy spectrum in turn has been shown to 
crucially influence the superconducting transition tem- 
perature and the critical current of the wire network. 
Before concluding this paper, we would like to make a 
few remarks. First, though we treated only the T3 ge- 
ometry in this paper, we expect that the field-induced 
localization properties of the T4 geometry 14 will undergo 
similar effects of field modulation. Second, for the sake 
of simplicity in the calculation, we dealt with a simple 
case where the period of field modulation is commensu- 
rate with the lattice period. The energy spectrum may 
exhibit more complicated band structures when the two 
length scales are incommensurate with each other. Fi- 
nally, note that the magnetic flux treated in this paper 
is neither uniform nor random but periodic. Thus, the 
effects of a random magnetic flux is still an open ques- 
tion, and it would be interesting to study whether or not 
AB cages remain bounded and/or how the localization 
properties of the system are influenced by introducing a 
random flux. Our naive expectation is that switching on 
a random flux at / = 1/2 will induce an energy band of 
localized states. 
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FIG. 1. A portion of the 2D bipartite periodic tiling. 
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FIG. 2. Energy spectrum as a function of / for T x = 3. (a) 
(3 = 0.0, (b) 13 = 0.3. 
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FIG. 3. Energy spectrum as a function of f3 for / = 1/2 
and T x = 3. 




FIG. 4. Density of states for / = 1/2 and T x = 3. (a) 
13 = 0.0, (b) 13 = 0.3, (c) f3 = 0.6, (d) (3 = 0.9. 
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FIG. 5. Energy spectrum as a function of (3 for / = 1/2 
and T x = 4. 
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FIG. 6. Density of states for / = 1/2 and = 4. (a) 
/3 = 0.0, (b) p = 0.3, (c) f3 = 0.6, (d) P = 0.9. 
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FIG. 7. Energy spectrum as a function of /? for / = 2/3 
and Ta, = 4. 




FIG. 8. 1 - T c (/)/T c (0) as a function of / for = 3. 
Continuous line: (3 = 0, dashed line: (3 = 0.3, dotted line: 
P = 0.6. 
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FIG. 9. l-T c (/)/T c (0) as a function of P for/ = 1/2. Con- 
tinuous line: T x = 2, long-dashed line: = 3, short-dashed 
line: T x = 4. 
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